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. Abstract. Recent attempts to fit Type la supernova data by modeling the 

Q I dark energy density as a truncated Taylor series have suggested the possibility of 

' metamorphosis, i.e., a rapidly evolving equation of state parameter, wde{z). However, 

we show that fits using that parametrization have significant problems. Evolution of 
wde{z) is both favoured and in some sense forced, and the equation of state parameter 
blows up or diverges in large regions of the parameter space. To further elucidate these 
\ problems we have simulated sets of supernova data in a A-universe to show that the 

' suggested "evidence" for metamorphosis is also common for wde = —1 

O 

o 
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1. Introduction 



^ \ Revealing the true nature of dark energy (DE) has become one of the most important 
?-j I tasks in cosmology. Considering the plethora of DE models proposed in the literature, 
^ \ a model independent reconstruction of DE would be an appealing alternative to testing 
all models separately. In two recent papers [H |2] attempts were made by Alam et al to 
^ ^ reconstruct the dark energy equation of state parameter wde{,z) in a model independent 
' manner, using the latest supernova data IH EI • In these two papers a truncated 
Taylor series was used to model the dark energy density Pde{z). The results indicate an 
evolution of wde{z), a behaviour they call metamorphosis. From the reported analysis, 
it would seem that this is a significant effect prompting for "exotic" models for the DE. 

Since other parameterizations |Sl El E! suggest that the Type la supernova data 
collected so far are consistent with the simplest DE model of all, that of a cosmological 
constant A {wde = —1), it is important to investigate how such different conclusions 
can be reached starting from the same sets of data. In this paper we argue that the 
method of model independent reconstruction proposed by Alam et al suffers from a 
number of serious shortcomings. For alternative methods of DE reconstruction see e.g. 
references [3 El CHI HH • 

A fundamental requirement on a model independent reconstruction of wde{z) must 
be that all DE alternatives are treated on an equal footing. Although the method of 
Alam et al at first sight seems to be capable of an exact reproduction of the equation of 
state parameter of the cosmological constant, it actually favours evolving wde{,z). The 
confidence contours, describing the level to which the reconstructed wde{.z) is known. 
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exhibit two related problems. First, contours enclosing regions of high confidence level 
(CL) and high redshift tend to diverge. Second, by construction, the low level CL regions 
shrink for high redshift s. 



2. Confronting supernova data with the model 

In this section we present the basic formulae needed for confronting the model of Alam 
et al with supernova data. The ansatz for the dark energy density proposed by Alam 
et al is a truncated Taylor series 

PDE = Aq + AiX + A2X^ (1) 

where x = 1 + z. In the following, we will assume that the universe is fiat and dominated 
by matter and DE. The Hubble parameter in such a universe, with the DE parametrized 
as in equation ((T)), is hence given by 

H{x) = HoiQMX^ + Ao + Aix + A-^x^fl^ . (2) 

If the matter density is assumed to be known, only two of the parameters describing 
the DE are independent. Ax and A2. The third parameter is given by the normalization 
condition, 1 = f2jv/ + + + y42. 

The equation of state parameter w{z) = p{z) / p{z) relates the pressure to the energy 
density and can be used to characterize an energy species. The dark energy equation of 
state parameter wde^z) can be calculated given the Hubble parameter and the matter 
density [12] 

, , (2x/3)d(\nH)/dx-l 

By using the parametrization in equation Q, wde{x) can be expressed in terms of the 
parameters Aq, Ai and A2 

^ ^ Ao + lAix + lA2X^ 

Using the derivative of equation ^ in the derivation of equation (j3)), is however a 
questionable step. If the parametrization of H{x) is incorrect, dH[x)/dx is likely to be 
even more erroneous. The parametrization is capable of reproducing the equation 
of state parameter for the cosmological constant {wde = —1, Ai = A2 = 0), and 
topological defects in the form of cosmic strings {wde = ^0 = ^1 = 0) and domain 
walls {wde = -|, Aq = A2 = 0). 

Models of DE are tested against observations by deriving distances to Type la 
supernovae through their measured brightness. Cosmic distances are related through 
the i^o^independent luminosity distance d'j^ to the brightness of a supernova measured 
in magnitudes 



m 



5\ogd'^ + M, (5) 
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where Ai = M + 25 — 51ogiJo will be treated as a "nuisance" parameter in the fitting 
procedure and the absolute magnitude of a supernova is denoted by M. The expression 
for the i^o^independent luminosity distance in a flat universe is 



Observed and theoretically predicted magnitudes can thus be compared by means of the 
equations (0) and © and the model of the DE enters via equation (j2)). However, since 
the integrand in (jH)) has to be real, i.e. VLmx'^ + + + A2x'^ > 0, only a part of 
the Ai — A2 parameter plane is allowed. As long as this condition is fulfilled, the dark 
energy density can be negative. 

3. Reconstructing the dark energy equation of state parameter 

In this section we reconstruct the dark energy equation of state parameter by the method 
proposed by Alam et al . 

We have used two sets of supernova data, which will be referred to as the small and 
the large data set. The small data set is a combination of supernovae from Tonry et al 
|H] and Barris et al |4 , and it is the same set as Alam et al used in reference The 
172 supernovae in the small set have low extinction Ay < 0.5 and z > 0.01. The large 
data set is an extension of the small set with 16 additional high redshift supernovae 
from Riess et al [S]. 

A maximum likelihood analysis of the data sets was used to estimate the best fit 
parameters of the dark energy model. The following negative log-likelihood function 
was used 



where rui and ai are observed magnitude and magnitude error of a supernova at redshift 
Zi- The parameter A4, related to the Hubble constant and the absolute supernova 
magnitude, was treated as a "nuisance" parameter in the fitting procedure described 
below. 

First, the negative log-likelihood function in equation was computed on a cubic 
lattice in the parameter space spanned by Qm, Ai and A2. Thereafter a Gaussian prior 

= 0.29 ±0.07, containing information on the matter density, was applied to £|. The 
Al — A2 plane containing the minimum of the negative log-likelihood function was then 
identified. Finally, the best fit values and the CL contours in this plane were translated 
into Wde{z) through equation (jH). 

The best fit parameters to the small data set {Ai = —6.68 and A2 = 2.84) and the 
99% and la CL contours in the parameter plane are shown in panel (a) in figure ^ The 
la CL contour is drawn in order to facilitate comparison with reference jl] and |2]. Note 
that Alam et al report their confidence regions as la levels, i.e. as Cmax ~ ^ in their 
I Note that the specific choice of prior of is not critical for our argument. 




(6) 
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maximum likelihood analysis. As two parameters are jointly fitted this corresponds only 
to a 39.3% CL. The reconstructed history of the dark energy equation of state parameter, 
shown in panel (b) in figure ^ is clearly consistent with a cosmological constant at the 
95%, but not at the 68% CL. There are however some peculiarities of the reconstructed 
wde{,z) that should be noted. The 95% and 99% confidence regions drawn in figure [T] 
diverge at large redshifts. This corresponds to the contour regions crossing the dark 
shaded region in the Ai — A2 plane, where wde{z) grows very rapidly. The 68% and la 
confidence regions, on the other hand, converge at large redshift. Moreover, there is a 
waist at z ~ 0.2 where all contours are squeezed together. These peculiarities indicate 
that this reconstruction is in fact not model independent, as will be further explained 
in the next section. 

Figure |2l shows wde{z) obtained from the large data set with 16 additional 
supernovae. This equation of state parameter is consistent with the cosmological 
constant even at the 68% CL. However, the confidence regions still exhibit the same 
behaviour as the ones obtained for the small data set. 

4. Analysis of the dark energy equation of state parameter 

The unusual behaviour of the results in figure ^ and |2 appear since the dark energy 
equation of state parametrization (jU has a few peculiarities, which we discuss in this 
section. 

4.1. Diverging equation of state parameter 

The dark energy equation of state parameter will obviously diverge whenever the 
denominator in equation (0)) vanishes. Regions in the parameter plane where W£,e{z) 
blows up > 15) or diverges somewhere in the redshift interval < ^ < 1.8 are 

shown in dark gray in e.g. panel (a) in figure Q The point Ai = A2 = 0, corresponding 
to the cosmological constant, lies very close to this region. If the best fit to data is near 
this point, the error in wde{z) will blow up. Diverging equation of state parameter is 
the reason why the 95% confidence region in e.g. figure ^b) blows up. 

4-2. High redshift limits and shrinking confidence regions 

From equation one can easily see that the parametrized equation of state depends 
on the order of the Taylor expansion at the high redshift limit. If, for instance, only 
the first order is considered, wde — | for high redshifts, while for the second order 
wde — |. Since the proposed ansatz ([T)) is a truncated Taylor series we cannot expect 
equation to be valid at very high redshifts. However, the effects of these high redshift 
limits are noticeable already at the moderate redshifts where we try to model DE. These 
limits are responsible for the shrinkage of the non-diverging confidence regions at high 
redshift (e.g., see the 68% confidence region in figure mb)). 
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Figure 1. The results of the reconstruction of DE from the small data set, consisting 
of supernovae from Tonry et al and Barris et al . Panel (a) shows the A\—Ai parameter 
plane, where the best fit value to the data {Ai = —6.68 and A2 = 2.84) is indicated by 
a star. The solid and dashed contours indicate the 99% and la confidence levels (CL) 
respectively. In panel (b) the reconstructed equation of state parameter corresponding 
to the best fit (solid line) and various confidence regions are drawn. The obtained 
equation of state parameter is consistent with a cosmological constant at the 95% CL, 
but not at the 68% CL. Note that Alam et al report their confidence regions as la 
levels, i.e. as Cmax ~ 5 their maximum likelihood analysis. As two parameters are 
jointly fitted this corresponds only to a 39.3 % CL. The shown discontinuities in the 
95% and 99% CL regions correspond to the contours crossing the dark shaded region 
in the Ai — A2 plane in panel (a), where wde{z) grows very rapidly. 
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Figure 2. The results of the reconstruction of DE from the large data set, consisting 
of the small data set and 16 additional high-z supernovae from Riess et al . In panel 
(a) the best fit values to this data set (^i = —1.44 and A2 = 0.61) is indicated by a 
star. The solid and dashed contours indicate the 99% and la confidence levels (CL) 
contours respectively. Both CL contours cross the region in the parameter plane where 
wde{z) blows up or diverges (dark gray). The equation of state obtained for this data 
set, shown in panel (b), is consistent with the cosmological constant not only at the 
95% CL, but also at the 68% CL. 



4-3. Stability of the dark energy equation of state parameter 

The dark energy equation of state parameter described by equation (jH) can be thought 
of as a family of curves, where each curve is parametrized by ^1 and A2. In this 
section we study the stabihty of these curves to perturbations in the parameters Ai 
and A2. Moreover, we will show that some curves in this family become more stable 
with redshift, with respect to these perturbations, while others become more unstable. 
The change 5wde in a curve, described by the parameters A\ and due to small 
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perturbations 6Ai and 6A2 is given by 

^DE - ^(A?, Al)5A, + ^(A?, Al)6A,. (8) 

The curve that corresponds to the cosmological constant {Ai = A2 = 0) becomes 
increasingly unstable due to perturbations as redshift increases, 

A curve that starts close to the line wde{x) = — 1 will thus departure from this line with 
increasing redshift. Equation (j3)) can also make an exact reproduction of two types of 
topological defects. The equation of state of domain walls {Aq = A2 = 0) grows stable 
to perturbations in Ai with z, but not to perturbations in A2 

~ on o ^ + on ^ ( + - ) "^^2- (10) 
3[1-Um>x 3[1-Um>\ xJ 

The curve corresponding to cosmic strings {Aq = Ai = 0) is stable to perturbations in 
both Ai and A2 at high redshifts 

The curves corresponding to topological defects are thus very special. All other curves 
(except the one corresponding to the cosmological constant) approach one of these 
two curves as redshift increases. The effect of perturbations in Ai and A2 depends 
on the value of these parameters. This is in contrast to perturbations of linear 
parameterizations of wde{z), e.g. wde{z) = Wq + Wiz, which are independent of 
parameter values. From the above discussion about the stabihty of wde{x) with redshift, 
we see that the parameterization in equation (jU treats the cosmological constant 
unfairly in favour of other DE alternatives. 



4.4- Forced evolution of the dark energy equation of state parameter 

The parametrized equation of state parameter will not only favour, but also "force" 
evolution. All curves, describing an evolving equation of state parameter, will cross at 
a certain redshift. At this redshift z^, = — 1, where x^, = —A1/2A2, the numerator 
and denominator are equal and the dark energy equation of state assumes the value 
u!de{z^) = —1. The parameters Ai and A2, obtained by fitting supernova data to the 
ansatz (^, are highly correlated and the ratio A1/A2 is thus nearly constant. The value 
of this ratio depends on the slope of the CL contours in the Ai — A2 plane which in 
turn is correlated with the assumed prior on the matter density Qm- This implies that 
the location of the point x^ = —A1/2A2 is more or less independent of the parameters 
Ai and A2 for a given Qm- The parametrization (jH) of the equation of state parameter 
thus favours DE evolution. If the reconstructed equation of state at present is less 
(larger) than minus one it must have been larger (less) than minus one in the past. 
Any significant departure from the line W£)e{x) = — 1 will hence force the reconstructed 
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equation of state parameter to cross that line. Since most curves will pass through this 
point, the family of curves will have a waist centered at = — 1/2A;, where k denotes 
the slope of the CL contours in the A1 — A2 plane. The transition redshift in e.g. figure^] 
is z^: ^ 0.2 which corresponds to A; = —0.42, in good agreement with the slopes of the 
ellipses in panel (a). 

5. Simulations 

To test the reliability of this way of reconstructing the dark energy we tried to recover 
a fiducial cosmology from simulated data sets generated using SNOC jlSj, a Monte- 
Carlo simulation package for high-z supernova observations §. We generated 500 data 
sets, each one consisting of 172 supernovae with the same distribution of redshifts 
and magnitude errors as the small data set used in section 01 A flat universe with a 
cosmological constant A and Qm = 0.3 was used as fiducial cosmology. The parameters 
Ai and A2 were fitted to each simulated data set with a fixed value of Qm- 

A scatter plot of best fit values to the simulated data sets is shown in figure El 
From this plot it is clear that the best fit values of the parameters Ai and A2 are highly 
correlated, as anticipated in section lOl Of the simulated data sets, 61.8% were found to 
be consistent with the cosmological constant at the 68% CL. These are indicated by full 
circles in the plot. Data sets not consistent with the cosmological constant at the 68% 
CL, on the other hand, are indicated by open circles. The best fit to real data obtained 
in section 121 is indicated with a star in the expanded figure in figure El The expansion 
shows a region where simulated data sets have a reconstructed dark energy equation of 
state parameter similar to the one obtained from real data. The simulated data sets in 
this region are inconsistent with the fiducial cosmology at the 68% CL. Figure El shows 
a few examples of wde{z) reconstructed from simulated data sets. The reconstructed 
wde{z) in panel (b) is almost indistinguishable from the one obtained from real data 
(see figure IH). Rapidly evolving equation of state parameters, resembhng the behaviour 
of the result obtained from real data, can thus be obtained from simulated data sets 
in a A cosmology {wde = —1) with the same distributions of redshifts and magnitude 
errors as the real set. 

We also simulated 500 data sets with the same distribution of redshifts and 
magnitude errors as the large data set used in section El Figure El shows a scatter 
plot of the best fit values to these simulations. In this case 64.9% of the simulated data 
sets were consistent with the fiducial cosmology on the 68% CL. As can be seen from 
figure El the 16 additional high redshift supernovae decrease the spread of the best fit 
values. 

To save computational time, a simplified fitting procedure was used in the 
reconstruction of wcie{z) from the simulated data sets, as Qm was kept fixed at the input 
value. The value of the matter density affects the slope of the confidence contour. By 
using a more realistic fitting procedure the spread of the best fit values would increase. 

§ The SNOC code is available at |ht t p : / / www . physto . se/ ~ ariel /snoc 
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Figure 3. A scatter plot of the best fit central values of Ai and A2 to 500 simulated 
data sets in a universe with a cosmological constant. The used cosmology corresponds 
to Ai = A2 = 0. The best fit values Ai and A2 are evidently highly correlated. 
The plot in the upper right-hand corner is an enlargement of the part inside the 
dashed rectangle. The filled circles represent simulated data sets consistent at the 68% 
confidence level with the fiducial cosmology. Data sets not consistent on this level 
are indicated by open circles. The star indicates the best fit value to real data. The 
dark energy equation of state reconstructed from data sets indicated by open circles in 
the neighborhood of the star, would thus look very similar to what we find from real 
data (see figure QJ. Examples of wde{z) reconstructed from a few of these data sets 
are shown in figure 0] 



6. Extended ansatz 

As we have seen above, equation (j3)) involves some difficulties associated with the 
reconstruction of the dark energy equation of state parameter. Could an extension 
of the ansatz overcome these difficulties? Studies of an arbitrary power series describing 
the dark energy density might help us to answer this question 

pDE = ^Akx\ (12) 

k 
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Figure 4. Examples of wde{z) reconstructed from three different simulated data sets 
with a cosmological constant are depicted in panel (b), (c) and (d). The distributions 
of redshifts and magnitude errors for these data sets were the same as for the small 
data set used in section |31 Best fit values to the three data sets are indicated in panel 
(a). The reconstructed equation of state parameter evolves rapidly in all three cases, 
which indicates that metamorphosis is common also in a A-universe. 



where k can be zero or any positive or negative integer. Perturbations in the parameters 
Ak will cause a change in p^)^; that is independent of these parameters 

W~5Za;'=Mfc. (13) 

k 

A description of the density of dark energy in terms of a power series will thus not favour 
any particular model. The equation of state generated by a power series is given by Q 

V A k 

If only one coefficient in p4|) is non-zero, wde{x) assumes a constant value. The power 
series parametrization of the equation of state parameter is hence capable of reproducing 




Figure 5. The same as in figure but the simulated data sets contain 16 additional 
supernovae at high redshifts. In this figure the star in the expanded figure indicates 
the best fit to the extended real data set. The additional supernovae have reduced the 
scatter of the best fit central values to the parameters Ai and A2 . 



a whole spectrum of constant wde{x) separated by Aw = | 

wde{x, ^ 0, Ak^j = 0) = -1 + ^. (15) 

The cosmological constant corresponds to j = 0. The high redshift hmit for a truncated 
power series with leading term resembles this spectrum 

Tl 

lim wde{x) = -! + -. (16) 

The change in wde due to any perturbation of the parameters A^ in the neighborhood 
of the cosmological constant {Aq 7^ and Akj^o = 0) is given by 

5wde ~ Y] ^x^5Ak. (17) 

The effect of perturbations in the parameters depends upon the leading term v4„x" in the 
power series. If n is positive, the change in the equation of state parameter will increase 
with redshift. However, if n is negative the change will decrease with redshift. In order 
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to discriminate between the cosmological constant and evolving DE alternatives, we 
have to contrast this behaviour with the curves describing evolution. The perturbed 
equation of state parameter corresponding to non-zero Ak at a high redshift (x 3> 1) 
with the leading term A^x"" is: 



The effects of perturbations decrease with redshift if n is positive, or increase with 
redshift if n is negative. This is the very opposite of the behaviour for the case of a 
cosmological constant. The instability problem discussed above can thus not be resolved 
by adding extra terms to the ansatz ((T)). 

7. Discussion 

The ansatz proposed by Alam et al may be useful for modelling the dark energy density, 
but its usefulness for revealing the nature of the DE seems limited. The parametrization 
of the dark energy equation of state parameter based on a truncated Taylor series 
involves a number of severe problems. Evolution is both favoured and forced by this 
parametrization. The cosmological constant is thus mistreated by the ansatz proposed 
in Ej. Not even an extension of the ansatz seems to be able to overcome these 
difficulties. 

The equation of state parameter expressed as in equation (jH) also diverges in large 
regions of the parameter space. More data, which focus the solutions to the stable 
regions, may solve this problem. However, the region close to the point describing a 
universe with a cosmological constant will be in a disfavoured part of the parameter 
space. 

The dark energy equation of state parameter reconstructed from the data sets 
presented in reference |2] and |1] is inconsistent with the cosmological constant at the 
68% confidence level. Simulations show that the rapidly changing behaviour of woEiz) 
could also be expected with this ansatz for a A universe. The scatter of the best fit 
parameters can be reduced if additional data points at high redshift are added to the 
simulated data sets. Our best fit to real data with 16 additional high redshift supernovae 
was consistent with the cosmological constant at the 68% confidence level. 

We concluded that the suggested "model independent" method of reconstructing 
the dark energy equation of state parameter in P, is in fact model dependent, and that 
the seemingly striking results are likely to be due to this deficiency. 
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